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Abstract 

The computation of the primordial power spectrum in multi-held inflation models requires 
us to correctly account for all relevant interactions between adiabatic and non-adiabatic 
modes around and after horizon crossing. One specific complication arises from derivative 
interactions induced by the curvilinear trajectory of the inflaton in a multi-dimensional field 
space. In this work we compute the power spectrum in general multi-field models and show 
that certain inflaton trajectories may lead to observationally significant imprints of 'heavy' 
physics in the primordial power spectrum if the inflaton trajectory turns, that is, traverses a 
bend, sufficiently fast (without interrupting slow roll), even in cases where the normal modes 
have masses approaching the cutoff of our theory. We emphasise that turning is defined with 
respect to the geodesies of the sigma model metric, irrespective of whether this is canonical 
or non-trivial. The imprints generically take the form of damped superimposed oscillations 
on the power spectrum. In the particular case of two-field models, if one of the fields is 
sufficiently massive compared to the scale of inflation, we are able to compute an effective 
low energy theory for the adiabatic mode encapsulating certain relevant operators of the 
full multi-field dynamics. As expected, a particular characteristic of this effective theory is 
a modified speed of sound for the adiabatic mode which is a functional of the background 
inflaton trajectory and the turns traversed during inflation. Hence in addition, we expect 
non-Gaussian signatures directly related to the features imprinted in the power spectrum. 



1 Present address: Theory Division, CERN, CH-1211 Geneva 23, Switzerland 



1 Introduction 



Single field slow-roll inflation [U [2] successfully accounts for many of the observed properties 
of the cosmic microwave background (CMB), including the near scale invariance of the power 
spectrum of the primordial density fluctuations that seed the observed CMB anisotropies [3]. 
Although one could claim that a large subset of the simplest models of single field inflation 
remain perfectly compatible with current CMB precision measurements [H [5], a direct and 
accurate reconstruction of the primordial spectrum from CMB data is still limited by various 
degeneracies in the priors and systematics adopted in our reconstructions [6-8J. It may certainly 
be the case that the CMB data implies the presence of various features in the primordial power 
spectrum other than the nearly scale invariant power law parametrization anticipated from the 
simplest single field slow roll models [SHE]- Upcoming data, such as that from the Planck 
satellite promises to provide new handles on the overall shape of the spectrum and, particularly 
in combination with other data sets, could help us determine the precise nature of any possible 
features in it. If present, such features could lead to quantitative new tests on the single field 
slow-roll paradigm |17[ [T8] and constitute strong evidence in favour of the existence of additional 
degrees of freedom present during the evolution of density perturbations as the universe inflated. 

One particularly compelling possibility which we wish to discuss in this report is of features 
in the spectrum generated by heavy (relative to the scale of inflation) degrees of freedom which 
do not necessarily decouple from the dynamics of the inflaton. Although the effects of massive 
degrees of freedom on the density perturbations are known to quickly dissipate during inflation, 
there are evidently still a number of contexts where features in the primordial spectrum due to 
heavy physics can survive. It is well understood, for example, that departures from a Bunch- 
Davies vacuum as the initial condition for the scalar fluctuations will result in oscillatory features 
in the power spectrum (see for example Refs. |19l [20]). Other contexts in which features are 
generated in the power spectrum involve particle production during brief intervals -much smaller 
than an e-fold- as the universe inflates. Examples of this include those situations where a massive 
field coupled to the inflaton suddenly becomes massless at a specific point in field space [2"T1 - 
[25] . Here it is the transfer of energy out of the inflaton field and the subsequent backscatter 
of its fluctuations off the condensate of created quanta that can result in features in the power 
spectrum, as well as in its higher moments |26[ [27] . Yet another context where such features 
have been shown to arise is in chain inflation, where instead of slowly rolling down a smooth 
continuous potential, the inflaton field gradually tunnels a succession of many vacua [28J. 

The purpose of the present report is to demonstrate the existence of, and to understand the 
general conditions under which features in the power spectrum result in the context of inflation 
embedded in a multi-scalar field theory (see Refs. [29H3T] for other recent discussions on this). 
For this we consider models of inflation where all of the scalar fields remain heavy except for one 
(the inflaton) which rolls slowly in some multi-dimensional potential. An effective field theory 
analysis tells us that in such scenarios, inflation should proceed in exactly the same way as in 
the single field case, with subleading corrections suppressed by the masses of the heavy scalar 
fields (see for example Ref. [32]). In this framework it is easy to take for granted that a simple 
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truncation of any available heavy degrees of freedom is the same as having integrated them out. 
However, it can certainly be the case that the adiabatic approximation is no longer valid at 
some point along the inflaton trajectory (e.g. due to a "sudden" turn that mixes heavy and 
light directions), and higher derivative operators in the effective theory are no longer negligible 
even as inflation continues uninterrupted. 

In various models of inflation in supergravity and string theory, the inflaton is embedded in a 
non-linear sigma model with typical field manifold curvatures of the string or Planck scale |33l - 
135] . In this type of scenario the inflaton traverses a curvilinear trajectory generating derivative 
interactions between the adiabatic and non- adiabatic modea^] |36H38| . In this context, it is 
straightforward to appreciate heuristically that a sudden enough turn can excite modes normal 
to the trajectory and non-trivially modify the evolution of the adiabatic mode. We will see that 
the net effect of this trajectory will translate into damped oscillatory features superimposed on 
the power spectrum - the transients after a sudden transfer of energy between the excited heavy 
modes and the much lighter inflaton mode, and the subsequent re-scattering of its perturbations 
off the condensate of heavy quanta that redshift in short time |^j 

A typical potential exhibiting such a curved trajectory is depicted in Figure [T] It can be 
appreciated that there is always a heavy direction transverse to the loci of minima determining 
the inflaton trajectory. We should emphasise however that the focus of this work is more general 
and that a curved trajectory in field space is not exclusively due to the shape of the potential, 
but also depends on the particular sigma model metric defining the scalar field manifold: on a 
particular curve the two can be transformed into each other by suitable field redefinitions. With 
this perspective, curved trajectories appear in any situation where a mismatch exists between 
the span of geodesies of the scalar field manifold and the actual inflationary trajectory enforced 
by the scalar potential through the equations of motion [JT] . The previously described situation 
is in fact generic of realisations of inflation in the context of string compactifications, where a 
large number of scalar fields are expected to remain massive but with their vacuum expectation 
values depending on the field value of the background inflaton [4"2"H5"T] . 

Limits of certain cases we wish to study in this report have been explored recently in seemingly 
different, but related contexts. In Refs. [52], [53] for instance, the effects on density perturbations 
due to a circular turn with constant curvature in field space was explored within a two-field 
model. There, it was concluded that such a turn could render non-Gaussian features in the 
bispectrum but would not generate features in the power spectrum. In another recent publi- 
cation [54] , the effects of a sigma model with non-canonical kinetic terms motivated by string 
theory was explored within inflationary models where one of the fields remained very massive. 
There, an effective theory was derived describing the multi-field dynamics, characterised by hav- 

2 Here, by adiabatic mode we refer to the mode which fluctuates along the inflationary trajectory whereas 
non-adiabatic modes correspond to those whose fluctuations remain orthogonal to the trajectory. We will also 
frequently denote them as curvature and isocurvature modes in this report. 

3 We also note the investigations of [391 !4Uj . where inflation in a putative string landscape is modelled using a 
random potential. Here, the background inflaton effectively executes a random walk, resulting in features at all 
scales in the power spectrum. 
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Figure 1: A generic example of a potential where turns happens while one of the fields remain much heavier 
than the other. 

ing a speed of sound for the fluctuations smaller than unity (and therefore indicating the possible 
departure from Gaussianity of the CMB temperature anisotropies) . In the framework we are 
about to discuss, both examples are just different faces of the same coin: while a non-canonical 
sigma model metric can always be made locally flat along a given trajectory generally generating 
contributions to the potential with a curved locus of minima, it is also possible to find a field 
redefinition which makes the loci of flat directions of the potential look straight at the cost of 
introducing a non-canonical metric. 

As we will shortly demonstrate, the parameter determining how relevant a local turn in the 
background inflaton trajectory is for the effective dynamics of the adiabatic mode is given by the 
departure from unity of the quantity e@ = 1 + 40q/(k 2 M 2 ), where <fio is the speed of the inflaton 
background field, k is the radius of curvature of the curve in field space and M is the mass of the 
direction normal to the trajectory. Keeping in mind that during slow-roll inflation, the inflaton 
velocity is given by 0o = V2eMp\H , with e being the usual slow roll parameter, it follows that 

= 1 + SeM^H 2 / (k 2 M 2 ) . Thus, even with M 2 » H 2 , if the radius of curvature describing 
the turn is small enough, significant imprints of heavy physics on the dynamics of the adiabatic 
mode can arise. More generally, whenever > 1, some amount of particle creation takes place 
that backreacts on the dynamics of the adiabatic mode. Let us not forget that in addition to the 
scale invariance of the power spectrum, single field slow-roll inflation predicts that the observed 
CMB temperature anisotropies seeded by the curvature perturbation satisfy Gaussian statistics 
to a high degree of accuracy [56]. Interestingly, in the class of models examined in this work, if 
the normal direction to the inflaton trajectory is sufficiently massive (M 2 > H 2 ), it is possible 
to compute an effective action for the adiabatic mode capturing the relevant operators of the 
full multi-field dynamics. This effective theory has the characteristic that the adiabatic mode 
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propagates with a speed of sound given by 



(1.1) 



and therefore becomes a functional of (the curvature of) the trajectory traversed by the in- 
flaton [5T]. Interestingly, this result has as a special case the particular context of Ref. |54j 
and indicates the presence of non-Gaussian signatures correlated with features in the power 
spectrum]^] 

We have organised this paper in the following way. In Section [2] we present the general setup 
and the notations used throughout this work. There, we will emphasise the need for using a geo- 
metric perspective to describe the evolution of the homogeneous background. Then in Section [3] 
we proceed to examine the perturbations of the fields around a time dependent background, and 
consider their quantisation and provide general formulas for the power spectrum. Our formalism 
allows us to consider situations beyond the regime of applicability of existing methods, such as 
trajectories with fast, sudden turns (regardless of whether the sigma model metric is canonical 
or non-canonical), and any other situations in which the masses in the orthogonal direction are 
changing relatively fast along the trajectory while still remaining much heavier than H 2 . In Sec- 
tion [4] we apply the previous results to the particular case of two-field models. We also deduce 
an effective theory valid in the case where the field normal to the trajectory remains heavy and 
compare the spectrum of this effective, single-field theory with those found in the full, multi-field 
computation of the spectrum. Then, in Section [5] we examine in detail the evolution of adiabatic 
and non-adiabatic modes along the curved inflaton trajectory and compute the power spectrum 
for various cases. Finally, in Section [6] we provide our concluding remarks. 



2 Basic considerations 



Let us start our study by recalling some of the basic aspects of multi-field inflation and by 
introducing the notations and conventions that will be used throughout this work. Our starting 
point is to assume the following effective four dimensional action consisting of gravity and a set 
of N scalar fields d> a : 



S 



-g d A x 



M, 



pi 



R 



1 



(2.1) 



Here R denotes the Ricci scalar constructed out of the spacetime metric g^ v with determinant 
g. Additionally, (j) a (a = 1, • • - A/") denotes a set of scalar fields spanning a scalar manifold M. 
of dimension Af, equipped with a scalar metric 7 a (,. The scalar fields may be thought of as 

4 In the course of preparing this manuscript, we note with interest the results of [55) . where their parametrization 
of the non-decoupling parameter of the isocurvature directions £ relates as a specific realisation of our analysis. 



This is easiest seen through comparing expressions (23) in Ref [53] with (4.271 or (C.9I here. We emphasise that 



the non-decoupling effect, ascribed in these two papers to the non-trivial sigma model metric, is simply due to 
the non-geodesic nature of the trajectory obtained by putting the heavy field at its local minimum. 
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coordinates on Ai with Christoffel symbols given by 

C = \l ad (d bldc + d clbd - d dlbc ) , (2.2) 

where d a are partial derivatives with respect to the scalar fields 4> a . In terms of these, the 
Riemann tensor associated with Ai is given by 

^ a bcd = d c Tl d - d d Tl c + T a ce T e db - T a de T e cb . (2.3) 

It is also possible to define the Ricci tensor as M. ab = M c a c& and the Ricci scalar M = -f ab R a b- We 
shall be careful to distinguish geometrical quantities related to the four dimensional spacetime 
and the AA-dimensional abstract manifold Ai. We should keep in mind that, typically, there will 
be an energy scale Aj^ associated to the curvature of Ai, and hence, fixing the typical mass 
scale of the Ricci scalar as R ~ ^m- manv concrete situations, such as the modular sector of 
string compactifications, the scale corresponds to the Planck mass Mpi. The equations of 
motion for the scalar fields are given by 

nr + Wd^d^ = v a , (2.4) 

where V a = j ab d b V. In what follows we discuss in detail the homogeneous solutions (j) a = (/>o(t) 
to these equations where the scalar fields depend only on time. The discussion in the sequel 
follows closely the analysis in [371 EH] ■ 



2.1 Background solution 

We look for background solutions by assuming that all the scalar fields are time dependent 
4> a = (f>o(t), and that spacetime consists of a flat Friedmann- Robertson- Walker (FRW) geometry 
of the form 

ds 2 = -dt 2 + a 2 (t)5 ij dx i dx j . (2.5) 
Later on we will also work in conformal time r, defined through the relation dt = a dr. In this 



background, the equation of motion (2.4) describing the evolution of the scalar fields is given by 

^ + 3H<p a + V a = 0, (2.6) 

where H = a /a is the Hubble parameter characterising the expansion rate of spatial slices, and 
where we have also introduced the convenient notation DX a = dX a + Ti X dffi. On the other 
hand, the Friedmann equation describing the evolution of the scale factor in terms of the scalar 
field energy density is given by 

where </>q = 7 a f></ ) o0o- Notice- that cj)o corresponds to the rate of change of the scalar field vacuum 
expectation value along the trajectory followed by the background fields. It is also convenient 
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to recall the following equation describing the variation of H 



which may be deduced by combining (2.6) and (2.7). By specifying the metric ^ a b and the 
scalar potential V, these equations can be solved to obtain the curved trajectory in A4 followed 
by the scalar fields. To discuss several features of this trajectory without explicitly solving 
the previous equations, it is useful to define unit vectors T a and N a distinguishing tangent and 
normal directions to the trajectory respectively, in such a way that T a N a = 0. These are defined 
as 



DT b DT C \ 1/2 DT a 



««=.»(*) ^TW-^-s-j -s-, (2.9) 

where sjv(t) = +lj denoting the orientation of N a with respect to the vector DT a /dt. That is, 
if sjv(i) = +1 then iV a is pointing in the same direction as DT a /dt, whereas if sjv(i) = — 1 then 
N a is pointing in the opposite direction. Due to the presence of the square root, it is clear that 
iV a is only well defined at intervals where DT a /dt ^ 0. However, since DT a /dt may become 
zero at finite values of t, we allow sj\r(i) to flip signs each time this happens in such a way that 
both iV a and DT a /dt remain a continuous function of t. This implies that the sign of sn may 
be chosen conventionally at some initial time U, but from then on it is subject to the equations 
of motion respected by the background^] In the particular case where M is two dimensional, 
the presence of sjv(i) in ( |2.9[ ) is sufficient for N a to have a fixed orientation with respect to T a 
(either left-handed or right-handed). This will become particularly useful when we examine two 
dimensional models in Section 01 

Observe that the tangent vector T a offers an alternative way of defining the total time deriva- 
tive D/dt along the trajectory followed by the scalar fields. This is: 

j t = ^ T a V a = ^oV^ . (2.10) 



Now, taking a total time derivative to T a , we may use the equation of motion (2.6) to write 

= -T Ta - T" ( 3 ^o + V a ) . (2.11) 

at (j) fa \ J 

Then, by projecting this equation along the two orthogonal directions T a and A^ a , we obtain 
the following two independent equations 

0o + 3^o + ^ = O, (2.12) 

^ = -SjV, (2.13) 
dt 0o 



5 We are assuming here that the background solutions 4> a = (f>o(t) are analytic functions of time, and therefore 
we disregard any situation where this procedure cannot be performed. 
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where we have defined = T a V a and Vm = N a V a to be the projections of V a = d a V along 
the tangent and normal directions respectively. It is not difficult to verify that V a lies entirely 
along a space spanned by T a and N a . That is, we are allowed to write V a = V^T a + V^Na. 
To anticipate the study of inflation within the present setup, it is useful to define the following 
dimensionless quantities: 

e = _ i L- ^ (2U ) 



1 Dd& 

V = = • (2.15) 

Hfo dt v ; 

We will not assume that these parameters are small until much later, where inflation is studied 



in the slow-roll regime (see Section 4.3). Similarly to the case of V a , the vector r\ a may be 



decomposed entirely in terms of T a and N a as 

7/ a =77||T a + 77 ± iV a , (2.16) 
- ^° (2.17) 



Vx , ( 2 .1 8) 

(bo H 



where we have used (2.6) to simplify a few expressions. Observe that r/j_ is directly related to 



the rate of change of the tangent unit vector T a , since (2.13) can be written as 

DT a 

-Hi 1± N a . (2.19) 



dt 



Comparison with (2.9) shows that sign(r/j_) = — sn- This is one of our main reasons for having 



introduced s^it) in (2.9): it allows us to keep r]± continuous and avoid some unnecessary 
difficulties encountered in the definition of isocurvature modes 

Moving on with this discussion, we can relate rj± to the radius of curvature k characterising 
the bending of the trajectory followed by the scalar fields. To do so, let us recall that given a 
curve 7(0q) in field space parameterised by d(f>o = <j>odt, we may define the radius of curvature 
At associated to that curve through the following relation: 

1 / DT b DT c ^ 1/2 



k \ dcpo dcpo ) 

Here k stands for the radius of curvature in the scalar manifold Ai spanned by the (j) a fields, 
and therefore it has dimension of mass. Figure [2] shows the relation between the pair of vectors 



6 In Ref. [30] for instance, a similar parameter rj± is introduced but with a fixed sign, which if the 'slow turn' 
approximation were to be violated (a regime we are mainly interested in), would produce a numerical overshoot 
in the evolution of curvature and isocurvature perturbations that has to be taken into account, and is naturally 
accounted for by our definition. 
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Figure 2: The figure shows schematically the relation between the tangent vector T a , the normal vector N a and 
the radius of curvature k. 



T a , N a and the radius of curvature k. Using (2.10) and comparing the last two equations we 
find the following relation between k and rj±: 



H\V±\ 



(2.21) 



By definition any geodesic curve j(<f>o) in A4 satisfies the relation Dcj) a /dt oc (f) a , which cor- 
responds to the case k~ 1 = 0, or alternatively, to the case i]± = 0. Thus, we see that the 
dimensionless parameter rj± is a useful quantity that parameterises the bending of the inflation- 
ary trajectory with respect to geodesies in M. It is interesting to rewrite the previous relation 
by replacing <j>Q = \^2eHMp\ coming from the definition of e presented in (2.15). One obtains: 



\V±\ 



Mpi 



(2.22) 



Then, if the radius of curvature is such that k Mpi, one already sees that rf^ S> 2e. We shall 
come back to this result later when we study curved trajectories in the slow-roll regime e <C 1. 
To continue, we may further characterise the variation of A^ a as: 



DN a 
dt 



1 



Hr]±T a + -—P ab V^V b , 



H V± 



(2.23) 



where we have defined the projector tensor P' 



ab 



7' 



ab 



T a T b — N a N b along the space orthogonal 



to the subspace spanned by the unit vectors T a and N a . That is, P ab N h = and P a bT b = 0. 
Details on how to obtain the previous relation can be found in Appendix |XJ Observe that in the 
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particular case where M. is two dimensional, one has 7 a b = T a Tf, + N a Nb and therefore P a b = 
identically. 

To finish this section, let us state some useful relations that will be used throughout the rest of 
this work. First, by using the definitions for e and 7711 in (2.15), we may rewrite the background 

V 



equations (2.7) and (2.12) respectively as: 

3 - 



With the help of (2.15) these two relations may be put together to yield: 

;uf>i / 1 ,, \ ~ 

e 



Next, by deriving (2.12) with respect to time and using the definitions for e and 



3(e + 77||) =M, 



1 , l -^(3-e) + ^||r ? ||, 



1 



(2.24) 
(2.25) 

(2.26) 
we deduce 
(2.27) 
(2.28) 



Both (2.26) and (2.27) are exact equations linking the evolution of background quantities with 
the scalar potential V. It may be already noticed that if e, 7711 and £|| are all much smaller 
than unity, then we obtain the usual relations for the slow roll parameters e and 7711 in terms of 
derivatives of the potential: 



e + 



v 



pr 



V 



(2.29) 
(2.30) 



We shall come back to these relations later, when we consider the evolution of the background 
in the slow roll regime. 



3 Perturbation theory 

The notation introduced in the previous section provides a useful tool to analyse perturbations 
5(j) a about the background solution <j) a = by decomposing them into parallel and normal 

components with respect to the inflaton trajectory. In what follows we proceed to study the 
evolution and quantisation of these perturbations. First, we consider scalar field perturbations 
by expanding about the background 4> a (t,x) = ^(t) + 5<p a (t,x). It is well known that the 
equations of motion for the perturbed fields can be cast entirely in terms of the gauge-invariant 
Mukhanov- Sasaki variables [571 EE] 
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where tp is the curvature perturbation of the spatial metric. The equations of motion for these 
fields are found to be [59] 

D 2 () a + 3R D<r _ Yl Q a + c a bQ b = ; (3 2) 



dt 2 dt a 2 



where V 2 = 5 l ^didj is the spatial Laplacian and where we have defined the tensor C a b as: 



C a b = V b V a - <j) 2 ] M. a cdb T c T d + 2eS- (T a V b + T b V a ) + 2e(3 - e)H 2 T a T b . (3.3) 

We notice here that C ab = j a cC c b is symmetric. It is convenient to rewrite the set of equations 
(3.2) in terms of perturbations orthogonal to each other. With this in mind, we introduce a 
complete set of vielbeins e J a = e^(t) and work with the following quantities: 

Q I (t,x) = e I a (t)Q a (t,x). (3.4) 

The a-index labels the abstract scalar manifold M whereas the /-index labels a local orthogonal 
frame moving along the inflationary trajectory. Recall that vielbeins are defined to satisfy the 
basic relations e I a e b 7 'y ab = S IJ and e^e^Sjj = ^ ab . From these relations one deduces the identities 

e iJ?- e « =- e «£. e i (3.5) 
a dt J J dt a y ' 

„ D j j D . . 

"a* < 3 - 6 > 

from which it is possible to read 

Q J =4^-y^Q J , (3.7) 

Q 1 =ei^ - 2Y I J Q J - (Y l K Y K j + Y 1 j) Q J , (3.8) 

where we have defined the antisymmetric matrix Yjj = —Yji as: 

T)p a 

= < 39 > 

Before writing down the equations of motion respected by the fields Q 1 , it is useful to notice 
that the matrix Yjj allows us to define a new covariant derivative T>/dt acting on quantities 
such as Q 1 labelled with the /-index in the following wa}Q 

V Q i = qI + y i jQ j (310) 
dt 



This definition allows us to rearrange (|3.7|) and (|3.8|) and simply write 

(3.11) 
(3.12) 



VQ 1 


iDQ a 


dt 


6a dt ' 


V 2 Q J 


I D 2 Q a 


dt 2 


~ 6a dt 2 



7 It may be noticed that we can write Y 1 j = (e^dbe^ + e^r^ c ej) </>q = uJb 1 j4v where uJb 1 j are the usual spin 
connections for non-coordinate basis, hence justifying the definition of the new covariant derivative of 13.101. 
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Thus, the equations of motion for the perturbations in the new basis become 



V2QI + 3H VQI 



V s 



dt 2 



dt 



Q 1 + C I J Q J = , 



(3.13) 



where Cjj = ei a e b jC a i,. To deal with the above set of equations, it is convenient to take one last 
step in simplifying them and rewrite them in terms of conformal time dr = dt/a, and a new set 
of perturbations v 1 = aQ 1 . These redefinitions induce a re-scaling of the covariant derivative 
(3.10) in the form D/dr = aD/dt, from where we are allowed to write: 



Vv 1 



dv 1 



+ Z*jv J 



(3.14) 



where Zjj = aYjj. With this notation at hand, the equations of motion for the -(/-perturbations 
are found to be 

v2 ' J vV + nV = o, (3.15) 



where O/j 



dr 2 

-a 2 H 2 {2 — e)6jj + a 2 Cu and we have used the definition of e to write a" /a 



a 2 H 2 {2 — e). For completeness, we notice that the previous equations of motion may be derived 
from the following action: 



S 



j drd 3 x Y, (jf) 2 ~ E( W ) 2 " ^ V 



(3.16) 



which can be alternatively deduced directly from the initial action (2.1 ) by considering all of the 
field redefinitions introduced in the present discussion. 



The set of equations (3.15) contains several non-trivial features. First, notice that the co- 



variant derivative T>/dr implies the existence of non-trivial couplings affecting the kinetic term 
of each field v . By the same token, under general circumstances the symmetric matrix 
does not remain diagonal at all times. In fact, as we are about to see in the next section, it is 
possible to choose to write this theory either in a frame where the M scalar fields are canonical 
(and therefore without non-trivial couplings in the kinetic term) , or either in a frame where fi/ j 
remains diagonal, but (in general) not both at the same time. 



3.1 Canonical frame 

Observe that by introducing the vielbeins e J a in the previous section, we have not specified any 
alignment of the moving frame. In fact, given an arbitrary frame, characterised by the set e^, it 
is always possible to find a canonical frame where the scalar field perturbations acquire canonical 
kinetic terms in the action. To find it, let us introduce a new set of fields u 1 defined out of the 
original fields v 1 in the following way: 

v 1 ^, x) = i? 7 j(r, n)u J (r, x) , (3.17) 
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where R 1 j(t, Tj) is an invertible matrix denned to satisfy the following first order differential 
equation: 

(3.18) 



d R 1 
d^ Rj 



-Z 1 rR K j 



with the boundary condition R 1 j(ji,Ti) = 5 1 j set at some given initial time Tj. Let us addition- 



ally define a new matrix S 1 j to be the inverse of R 1 j, i.e. 



S i kR k j 



R i kS k ^ 



it is possible to see that S 1 j satisfies the following similar equation 

d 



dr 



Sf 



-Z kSi 



K 



S 1 j. Then, 



(3.19) 



where we used the fact that Zjj = —Zjj. Since both solutions to (3.18) and (3.19) are unique, 
then the previous equation tells us that Sjj = Rji, that is, S corresponds to R T the transpose 
of R. This means that for a fixed time r, Rij(t, t{) is an element of the orthogonal group O(AT), 
the group of matrices R satisfying RR T 
symbolically written as 



1. The solution to (3.18) is well known, and may be 



R(T,r i ) = t + J2 



(-1)' 



n=l 



nl 



T[Z(n)---Z(T n )]d n T = Texp 



drZM 



(3.20) 



where T stands for the usual time ordering symbol, that is T\Z{t\)Z{t2) ■ ■ ■ Z(r n )] corresponds 
to the product of n matrices Z(n) for which r\ > T2 > • • • > r n . Coming back to the u^-fields, 
it is possible to see now that, by virtue of (3.18) one has: 



dr 
dr 2 



-R'j 



dr 
d 2 u J 
dr 2 



(3.21) 
(3.22) 



Inserting these relations back into the equation of motion (3.15) we obtain the following equation 
of motion for the u^-fields: 

d 2 u ! 



dr 2 



VV + [R T {t) n R(t)] 1 ju J = . 



Additionally, it is possible to show that the action ( |3.16| ) is now given by 



S 



drd 3 x 



i.i 



u'u 3 



(3.23) 



(3.24) 



Thus, we see that the fields u 1 correspond to the canonical fields in the usual sense. This result 
shows, just as we have stated, that it is always possible to find a frame where the perturbations 
become canonical, but at the cost of having a mass matrix [R t (t)SIR{t)] ij with non-diagonal 
entries which are changing continuously in time. Another way to put it is that, while both 
R t (t)Q,R(t) and Q share the same eigenvalues, as long as R(t) varies in time, their associ- 
ated eigenvectors will not remain aligned. To finish, let us notice that by construction, at the 
initial time Tj, the canonical fields u 1 and the original fields v 1 coincide n^(rj) = v 1 ^). How- 
ever, it is always possible to redefine a new set of canonical fields by performing an orthogonal 
transformation of the fields. 
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3.2 Quantisation and initial conditions 



Having the canonical frame at hand, we may now quantise the system in the standard way. 
Starting from the action (3.24) it is possible to see that the canonical coordinate fields are given 
by u 1 whereas the canonical momentum is given by 11^ = du 1 J dr. To quantise the system, we 
demand this pair to satisfy the commutation relation 

[u\r t x),Hi(T,y)] = i5 IJ 6®(x-y), (3.25) 



otherwise zero. With the help of the R transformation introduced in (3.17) we can rewrite this 
commutation relation to be valid in an arbitrary moving frame. More precisely, we observe here 
that we are allowed to define a new pair of fields v 1 and 11^ given by 

v 1 =R I JU J 

V T 



(3.26) 
(3.27) 



n£=— v 1 =R I j(r,r i )Ui. 

From (3.25), this new pair is found to satisfy similar commutation relations: 

[« J (r,aj),I^(r,y)] = id IJ 6^(x - y) . (3.28) 

Following the convention, it is now possible to obtain an explicit expression for v J (x, r) in terms 
of creation and annihilation operators^] For this, let us write u J (x, r) as a sum of Fourier modes: 

d 3 k 



V ( T, X ) 



ikx 



t/(fc,r) 



3 *- £ [vi(k, r) a a (k) + vi*(k, r) a^-k) 



(3.29) 



(2tt) 3 / 2 
d 3 k 

(2^) 3 / 2 ' 

In writing the previous expression we have anticipated the need of expressing the fields ?/(t, x) 
as a linear combination of AT time-independent creation and annihilation operators al(fc) and 
a a (k) respectively, with a = 1, • • • N. These operators are required to satisfy the usual relations 

5 a ^ 3 \k-q), 



a a (k),al(q) 



(3.30) 



otherwise zero. This set of operators defines the vacuum |0) of the theory by their action 
a a (fe)|0) = 0. Since the operators a/ a (k) and a a (k), for different values of a, are taken to be 
linearly independent, then the time-dependent coefficients v^(k,T) appearing in front of them 
in (3.29) must satisfy the equation of motionj^] 

V 2 

7 v I a {k,r) + k 2 v I a {k,T) + n I jv J a (k,T) = 0. (3.31) 



dr 2 



8 From this point on, we continue working with the more general u I -fields instead of the canonical w 7 -fields. 
Nevertheless, we emphasise that the it^-fields allowed us to find the correct quantisation prescription for the 
i/-fields, which otherwise would not have been properly justified. 

9 It is crucial to appreciate that the Greek indices a label scalar quantum modes and not directions in field 
space, as capital Latin indices do. Different a-modes may contribute to the same fluctuation along a given 
direction /. The quantities linking these two different abstract spaces are the mode functions v a (k, r) whose time 
evolution is dictated by ( |3.3f \ . A similar scheme to quantise a coupled multi-scalar field system may be found in 
Ref. [50]. 
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Observe that there must exist Af independent solutions v^(k, r) to this equation (see Appendix |b| 
for a detailed discussion on the v a (h, r)-functions). 

Of course, a critical issue here is to set the correct initial conditions for the mode amplitudes 



v^(k, t) in such a way that the commutation relations (3.28) are respected at all times. As a first 
step towards determining these initial conditions we notice that at a given initial time r = Tj 
we may choose each mode v^(k,T) to satisfy the following general initial conditions: 



vl{k,Ti) = e I a v a (k) , 
1 Vn, -(k,Ti) =e I a -K a {k) 



dt 



where is a complete set of unit vectors satisfying Sije^e^ 



la? 



and 5 a/3 e I a en 



(3.32) 
(3.33) 

5 IJ , which 



should not be confused with the vielbeins defined in (3.4), and v a (k) and n a (k) are factors 
defining the amplitude of the initial conditions. In order for the commutation relations to be 
fulfilled, these initial conditions must satisfy: 



v a (k)-n* a (k) - v* a {k)-n a (k) = i 



(3.34) 



which are the analogous relations to the Wronskian condition in single field slow-roll inflation. 
Since the operator T>/dr = d/dr + Z mixes different directions in the ■iZ-field space and since 
in general the time-dependent matrix Qjj is non-diagonal, then the mode solutions v^(k,r) 



satisfying the initial conditions (3.34) will not remain pointing in the same direction (nor will they 
remain orthogonal) at an arbitrary time r ^ Tj. In Appendix |B] we show that the commutation 



relations of (3.28 ) are consistent with the evolution of the v^(t, k) dictated by the set of equations 



of motion (3.31 ) 



In the previous expressions the set of unit vectors e„ are arbitrary. Moreover, the amplitudes 



v a (k) and n a (k) entering (3.34) are in general not uniquely determined, as there is a family 
of solutions parameterised by the relative phase between v a (k) and ir a (k). Indeed, notice that 
without loss of generality we may write 



ir a (k) 



,(*) 



2v* a (k)sm6 a (k) 



(3.35) 



where 6 a {k) is a set of real phases relating both amplitudes. Any value for 6 a (k) will satisfy the 



commutation relations (3.28), and therefore they specify different choices for the vacuum state 
|0). Although in general it is not possible to decide among all the possible values for 8 a (k), 
fortunately, in the context of inflationary backgrounds a — > as r — > — oo and a particular 
choice for these phases becomes handy. Indeed, observe that in the formal limit a — > one has 
Z[j — > and tlu — > 0, which is made explicit by (3.14) and (3.15), and the equations of motion 



(3.31 ) become: 



d^ 



+ k 2 )v I a (k,T)=0. 



(3.36) 
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In this limit there is no mixing between different a-modes and perturbations evolve as if they 
were in Minkowski background^] In this case, we are free to choose = 5^ and the solutions 



to (3.36) satisfying the commutation relations (3.28) may be chosen as: 



1 



-ikr 



2 k 



(3.37) 



Thus we see that in the limit a — > (r — > — oo) we may choose modes in the Bunch-Davies 
= it/2. We will come back to these conditions in Section [5] where we set initial 



vacuum 



conditions on a finite initial time surface where (3.37) cannot be exactly imposed. 



3.3 Two-point correlation function 

To finish this general discussion on multi-field perturbations, we proceed to define the spectrum 
for the perturbations v 1 (t, x). The power spectrum, the Fourier transform of the two-point 
correlation function, is defined in terms of the Fourier modes as 

<0 |«'(fe,T)t; J *(g,r)| 0) = 6®(k - q^V^ (k,r) . (3.38) 
In terms of the mode amplitudes v^(t, k), this is found to be 

Vv'ik, r) = £ vi( T , k)v J a * (k, r) . (3.39) 

a 

Since the commutation relations require J2 a [v^(k,T)v^*(k,r) — v^(k, r)v^(k, r)] = (see Ap- 
pendix [bJ we see that the spectrum Vl J is real, as it should be. Additionally, the two point 
correlation functions in coordinate space may be computed out of Vy J as: 

<0|^(r,^ J (r,y)|0> = — J —Vi J (k,T)e- lk < x -yK (3.40) 

We may also define the power spectrum associated to the Q fields instead of the v fields. 
Recalling that Q 1 = v 1 /a, the power spectrum for these fields at a given time r is then given 
by: 

a 

This expression will be used to compute the power spectrum of the curvature perturbation 
produced during inflation. Although, in this section we have chosen to exploit a notation whereby 
Greek indices a label quantum modes, notice that this formalism is equivalent to the use of 
stochastic Gaussian variables, as in Ref. [61] (see also Ref. |49|). 



3 To be more rigorous, in inflationary backgrounds this limit is obtained for fc-modes such that their wavelength 



is much smaller than the de Sitter scale k 2 ^> a 2 H 2 . 
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4 Models with two scalar fields 



We now study the evolution of perturbations in systems containing only two relevant scalar 
fields. In this case, it is always possible to take the set of vielbeins {ej} to consist entirely in 
= T a and = N a defined in Section 2.1 Then, the projection tensor P a b introduced in 



(2.23) vanishes identically and one is left with the following relations: 

DT a 



Hrj±N a , (4.1) 



dt 
DN a 

— =H V± T a . (4.2) 
At this point we notice that the normal vector N a has always the same orientation with respect 



to the curved trajectory, which is due to the presence of the signature function sn in (2.9). For 
definiteness, let us convene that the normal direction N a has a right-handed orientation with 
respect to T a as shown in Figure [3} With this convention rj± changes signs smoothly in such a 
way that if the turn is towards the left then r]± is negative, whereas if the turn is towards the 
right then r/± is positive. A concrete choice for T a and N a with these properties are: 

' ' - (4.3) 



00 

N a = -J— (-722^ - 7120 1 , 7H0 1 + 7210 2 ) , (4.4) 
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where 7 = 711722 — 712721 is the determinant of 7^. To continue, parallel and normal pertur- 
bations with respect to the inflationary trajectory are then given by 



v 

,.A f 



aQ T = aT a Q a 



aQ N = aN a Q a 



(4.5) 
(4.6) 



By choosing this frame, one finds that Ztn = —Znt = aHr]±. The coupled equations of motion 
describing the evolution of both modes v^(k, r) and (k, r) become: 



dr 2 
d 2 v% 
dr 2 

where we have defined: 



d 2 v T dv N 
9L _i_ nr a 

' ' " S dr 

dv? 



ev T a + ^ + n TN v^ + {n TT + k")vl 



C 2 v N 



dr 
dr 



0, 
0, 



C = Z TN = aHr]± . 



(4.7) 
(4.8) 

(4.9) 



In the previous equations, the symmetric matrix Qjj is defined in (3.3) and (3.15) and consists 
of the following elements: 



n TT = - a 2 H 2 (2 + 2e - 3r?, , + rj^ - 4e V]] + 2e 2 - rfi) , 
n NN = - a 2 H 2 (2 - e) + a 2 M 2 , 
n TN =a 2 H 2 V± (3 + e - 2r)\\ - £]_) , 



(4.10) 
(4.11) 
(4.12) 



where M = Vatat + J? Mpj e E is the effective squared mass of the t> -mode and R = 2Mj'jv'TA r = 

T a N b T c N d R 

abed is the Ricci scalar parametrsing the geometry of A4. Additionally, we have 

defined 



(4.13) 



Details on how to arrive at this specific form of Cljj for the case of two-field models are given in 
Appendix [Aj 



4.1 Power spectrum 



Expressions (4.7) and (4.8) consist of the equations of motion necessary to deduce the generation 



of the curvature perturbation in the case of two-field inflation. Once the solutions of the fields 



aQ 1 and v = aQ N are known, it is possible to define the curvature and isocurvature 



perturbations as 



9o 

S^Q N 

9o 



(4.14) 
(4.15) 
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respectively. Using equation (3.41) with I = J = T, the resulting power spectrum for adiabatic 
modes are found to be 

Vn(k,r) = ^V% T (k,T) = k l Y, v^(k,T)vl*(k,r), (4.16) 

4>0 47F a IV1 Pl e Q= l,2 

where a and e = (j) 2 / '(2Mp,iI 2 ) are functions of r. We can also compute the power spectrum for 
isocurvature modes and cross correlation as [361 E2] 



V s (k,r)=S-r^ N (k,r), 
VKs(k,T)=~V% N (k,T), 



(4.17) 
(4.18) 



respectively. They can give rise to observable signatures in the CMB power spectrum [62], but 
it depends on post-inflationary processes thus we do not consider them here. In this work we are 
primarily concerned with the computation of the power spectrum of the curvature perturbation 
1Z at the end of inflation. This corresponds to the quantity 



Vn{k) = T>n(k,T end ) , 



(4.19) 



where r en d is the time at which inflation effectively ends 11 But the computations of Vs and 
Vtis can be done in an identical way. 



4.2 Effective Theory 

If a hierarchy of scales is present in the matrix f2/j, then we can compute a fairly reliable effective 



theory out of the system (4.7) and (4.8). Indeed, by assuming that Qnn remains positive at all 
times and that 

V NN >|^tt| , (4.20) 

n NN ^>\n TN \ , (4.2i) 

then we may integrate the heavy mode v N out of the system of equations. By examining the 
specific shape of the entries 0,tt, &nn and FLtn we see that a generic requisite for this hierarchy 
to exist is 

M 2 » H 2 , (4.22) 
where M 2 is the effective mass of the heavy mode v N given by: 

M 2 = V NN + H 2 M^eR. (4.23) 

11 Since in multi-field inflation the adiabatic mode 1Z (as well as other background quantities) may continue 
evolving on super horizon scales, here we do not follow the standard practice of evaluating the power spectrum 
at horizon crossing time k = aH [64]. See Ref. [65] for a discussion of this point. 
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To compute the effective theory we proceed in the same way as in Ref. |41j . We focus on the 
mode a associated to slower oscillations due to the hierarchy. Omitting the a label, this mode 
is necessarily such that 

d 2 v N 



dr 2 



< a 2 M 2 v N 



(4.24) 



This allows us to disregard the second derivative of v in (4.8), and write v N in terms of v T : 

1 



n 



NN 



d( T 



(4.25) 



This expression for v N can be inserted back into the remaining equation of motion (4.7) to 
obtain an effective equation of motion for the light adiabatic mode v T . Then, by defining a new 
field ip as 



Hi + 47/j 



2 + e- 



vl + 



A: 2 



(Ml 

we finally arrive at the following effective equation of motion^ 2 ] 

ip" + e-^ T > k2) k 2 v + n(r, k 2 )<p = 
where the time dependent function Q(t, k 2 ) is found to be: 



n(r,k 2 ) =n (r) 



2 



aH^(l + €-77|i), 



Q (t) =-a 2 H 2 (2 + 2e 



(4.26) 
(4.27) 

(4.28) 

(4.29) 
(4.30) 



3??|| -4er?|| +Z\\T}\\ +2e 2 ). 

Notice that Qq, is precisely the mass term appearing in the conventional equation of motion 
for adiabatic fluctuations in single field slow-roll inflation. Furthermore, we note that in the 
case where the mass M approaches the cutoff of our theory, our results can be derived from an 
effective action for the adiabatic mode given by the action 

2 



drd 3 x 



dip 
drr 



V(p e 



-/3(r-V 2 ) 



Vip-<p Sl(r, -V 2 )<£ 



(4.31) 



where /3(r, -V 2 ) and fi(r, -V 2 ) are the functions defined in M.26h and (^4.29^ but with k 
replaced by —V 2 . This result corresponds to the generalisation of our previous work [5T] to the 
case of a slowly rolling background in the presence of gravity. A slightly more formal deduction 
of this effective theory may be found in appendix [Cj where we see that it can be viewed as a 
leading order effect at the loop level, and as such contains the higher dimensional corrections 
implied by the general arguments made in Refs. |32|, I66J . In Section [5] we shall compare the power 
spectrum obtained using this effective theory with the one obtained from the full set of equations 
for the perturbations. We anticipate that this effective theory is very reliable regardless of how 
large the values of f3 are. 



12 Please notice that eqs. (4.281 and (4.301 are corrected versions of those appearing in the published version of 
this article. 
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4.3 Slow-roll inflation in two-field models 



So far we have not assumed the slow evolution of background quantities. We now proceed to 
discuss the case of inflation realised in the slow-roll regime, where the scale of inflation H varies 
slowly. Our main interest is to study the effects appearing from curved inflationary trajectories, 
where r/± is non-vanishing. We will assume that the radius of curvature k may take values 
smaller than Mp\, corresponding to turns of the trajectory taking place at field scales smaller 
than the Planck scale. This situation is certainly allowed and depending on the value of e, it 
may render large values of rj± (recall (2.22) relating rj± and n). By the same token, we will 
consider models where the normal mode v has a large effective mass M 2 3> H 2 . 



4.3.1 Slow-roll parameters 



In general, given the background equations of motion (2.6), (2.7) and (2.8), we say that a given 



background quantity A is slowly rolling if its variation satisfies 

1 dA 



\Sa\ 



HA dt 



< 1. 



(4.32) 



5^, and therefore both H and 0o evolve slowly 
2 l(n »>r2 jy2^ 



Observe that we can write e = §h and r/|| 

if e <C 1 and \rju\ <C 1 respectively. Since e = cp 2 / (2Mp l H z ) then the condition |r?||| <C 1 also 
guaranties that e will remain varying slowly during inflation. It is useful to introduce a single 
small dimensionless number 8 <C 1 parametrising the slow-roll expansion 13 and demand that 
any quantity A to which slow-roll is imposed, generically satisfies 

1 dA 



0(5) 



(4.33) 



HA dt 

0{8). In the absence of clear evidence of it, for simplicity we 



which means e = 0(5) and rj 
shall not consider here hierarchies between different slow-roll parameters. Recall that (2.26) and 



(2.27) are exact equations relating the parameters e, 7711 and £|| to the shape of the potential V 



along the inflationary trajectory. Now, provided that all of these parameters are small, we may 
re-express these equations to leading order in 5: 



+ e =M, 



pr 



V 

Y± 
v 



(4.34) 
(4.35) 



These are the usual equations defining the slow-roll parameters in terms of the shape of the 
first and second derivatives of Vj^] As long as e <C 1 and \r)\\\ 1, the background geometry 

13 Current observations indicate that the order of such a reference parameter is given by the departure of the 
spectral index from unity S ~ \n-iz — 1|. We remark here for completeness that the e and rj parameters above 
correspond to the Hamilton- Jacobi slow roll parameters. 

14 Let us recall that the parameter 77 was originally introduced in the study of single field slow-roll inflation |67j 
as rj = Ml x V" IV. Therefore, in order to compare the present results with those following the original convention, 
we must write f] = 7711 + e. 
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evolves slowly and the scalar field velocity is determined by the attractor equation of motion 
3H(/)q + Vj, = 0. For completeness, notice from the definition of r/± in (2.18) that it is possible 
to write rj± = Vn / ' {y/2eM-p\H 2 ) . Then, using (4.35) we deduce 



2 n I Vn 



(4.36) 



which is valid to leading order in 5. This equation nicely relates the slope of the potential 
along the tangential direction T a with its counterpart Vn along the normal direction iV a . 



4.3.2 Perpendicular dynamics 



Let us now turn our attention to parameter rj± defined in (2.18). Notice that this parameter is 



not related to the slow-roll variation of any given background quantity A in the sense of (4.32) 



and therefore it is not constrained to be of 0{5). Moreover, (2.22) tells us that r/j_ may be 
large compared to 5 provided that the radius of curvature k is small compared to \^2eMp\. It 
is important to recognise that the curved inflationary trajectory ^ 0) has its origin in 

both the shape of the scalar potential V and the geometry of the scalar manifold where the 
theory lives. In particular, since H and <f>o are assumed to evolve slowly, then we expect the flat 
inflationary trajectory to remain close to the locus of points minimising the heaviest direction 
iV a of the potential. In other words, to ensure a bending of the trajectory we consider models 
where the potential is such that 

V NN » |V^| . (4.37) 

It is entirely clear that in the event that the inflationary trajectory is suffering a turn, it will 
not coincide exactly with curve minimising the heaviest direction, which is made explicit by the 
result Vn 

0, with M 2 given by 



real minima Vn\ 



r)±(j)oH found in (2.18). It is in fact easy to show that the departure A from the 
= is roughly given by the condition Vn + M 2 A 



(4.23). Then, with the help of (2.18) one finds that the ratio between the deviation A and the 



radius of curvature k is given by 



A 2 H 2 
— - V± 

K 



M 2 " 



(4.38) 



Observe that A/k is essentially the combination e 



°—l defined in (4.26) in the regime k 2 <C a?H 2 . 

is giving us 



Thus the parameter (3 appearing in the effective theory deduced in Section 4.2 
information regarding the dynamics perpendicular to the inflaton trajectory. 

It is important to check out whether the bending interferes with the flatness of the potential 



as felt by the adiabatic mode v T . Observe from (4.7) and (4.10) that the effective mass m 2 (r) 
of v T is given by 



m 2 (r) = Q T t ~ C 



-a 2 H 2 (2 + 2e 



3»7||) , 



(4.39) 



where we have neglected terms of 0(5 2 ). Note that m 2 (r) = ^o( r )> where ^o( r ) is the effective 
mass encountered in the effective theory deduced in Section |4.2| Thus, we see that ij± does not 
directly spoil the flatness of the potential V. Of course, one should explicitly verify in which 
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way a bending affects the value of e and r/|| by examining the evolution of the background. We 
however point out that there is no reason a priori for which fast and sudden turns with large 
values of rj± are not possible while staying in the slow-roll regime. 



4.3.3 Equations of motion in the slow-roll regime 



Putting all of the previous results together back into the set of equations (4.7) and (4.8), 
and neglecting terms of 0(5 2 ), we finally arrive at the following equations of motion for the 
perturbations and '■ 

d 2 vl . „ dv" 



dr 2 



2aHrj A 



dT 



+ a 2 H 2 



i 2 H 2 



2 - 2e + fy^ x£ + 2a 2 H 2 7] ± (2 - £j 



2„JV 



dT 2 



2aHr}_ 



■ dr 



+ a 2 W 



k 2 M 2 
a 2 H 2 + H 2 



2 + e 



,N 



+ 2a 2 H 2 r) ± vl 



(4.40) 
(4.41) 



where £j_ was defined in (4.13). In the next section we deal with these equations numerically 



for suitable choices of the background parameters, and compare the obtained power spectrum 



with that of the effective theory obtained in Section 4.2 We shall see how features in the power 
spectrum appear as a consequence of curved inflationary trajectory. 



5 Features in the power spectrum 

We now study the evolution of perturbations and analyse how features in the primordial spec- 



trum are generated along curved trajectories. To this extent, we solve (4.40) and (4.41) numer- 



ically for different background solutions representing curved trajectories and obtain the mode 



solutions which, with the help of (4.19), provide us the desired power spectrum at the end 
of inflation. For definiteness, we consider models of inflation with an inflationary period of at 
least 60 e-folds and set the initial conditions a few e-folds before this period starts. To avoid 
unnecessary complications with initial conditions, we considered models where turns in the tra- 
jectory only happen within the last 60 e-folds. Before this period, rj± = and the equations of 
motion determining the evolution of perturbations reduce to 



d 2 vl 


+ a 2 H 2 ( 


f k 2 


dT 2 


K a 2 H 2 


d 2 v% 


+ a 2 H 2 


( k 2 


dT 2 


\a 2 H 2 



+ 



2-2e + 3r/|| 
M 2 



H 2 



2 + e 



.N 



=0. 



(5.1) 
(5.2) 



Then, as long as e and ?7|| are small, we are allowed to make use of initial conditions (3.32 ), (3.33 ) 
and (3.34) with e„ = <5^, and v\(k) and V2{k) given by 



vi(k) 
v 2 (k) 



7T 



V(l - e)aiH. t 



. e *i("i+i)#(i) 



k 



7T 



=^1(^ + ^(1) 

V(l - e)<nHi U2 \(l-e)oiHi 



(1 - e)a,iHi 
k 



(5.3) 
(5.4) 



22 



where H„(x) denotes the first kind Hankel function, whereas Gtj and Hi are the values for the 
scale factor and Hubble parameter at the initial time Tj. Similarly, the quantities irx(k) and 



7T2(/c) entering the initial conditions (3.34) are given by the time derivatives of the previous 
expressions. On the other hand, the parameters v\ and vi are respectively given by: 



(3 



4(1 



3(r? 



(1 



(3 



M 2 



4(1 



(5.5) 
(5.6) 



Note that in the short wavelength limit k 3> a?,-ffi, the previous conditions match the mode 



fluctuations about a Bunch-Davies vacuum (3.37) discussed in Section 3.2 In all of the cases 



examined, we consider inflationary trajectories where e, r/n and £|| remain small during the 
interval of interest, while allowing different types of time variation of rj±, which is the quantity 
that parameterises the bending. 



5.1 Constant radius of curvature 

Let us start by considering the simple case in which r]± is constant during the whole period 
of inflation where currently accessible modes were generated. As we have already emphasised, 
if e remains nearly constant, then a constant tj± corresponds to a trajectory with a constant 
radius of curvature k. We find that the overall effect of having a constant turn is simply to 
normalise the amplitude of the spectrum, without modifying the usual single field dependence 
of the spectral index n-n in terms of the slow-roll parameters e and r]\\ (see also [521153]): 

n^-l = 2r ? ||-4e. (5.7) 

In the case M 2 /H 2 ^> 1, the predicted power spectrum obtained by the effective theory is 
indistinguishable from the one obtained by solving the full set of equations. Moreover, with the 
help of this effective theory, it is in fact possible to infer a simple relation between the power 
spectrum Vn{k) with r]± / and the analytical power spectrum V^\k) computed with r]± = 0. 
To this extent, notice that although f3(k, r) is a function of k, we see that after the physical 
wavelength of the mode becomes larger than the scale M _1 (i.e. k 2 /a 2 < M 2 ), the parameter 
(3(r,k) becomes effectively k independent, and we can write 

e P = l + W±^- (5-8) 

Since M 2 ^> H 2 , this happens before horizon crossing, and the relevant dynamics is well de- 
scribed by this ^-independent form of /3. Then the relation between Vn(k) and V^\k), as 
predicted by the effective theory, becomes: 

Vn{k)= (l + ^i^)rS\k). (5.9) 
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This result modifies the usual normalisation condition of the spectrum coming from the COBE 
data, leading to the following relation among the various parameters: 

(l + 4? ?i^) ?4 Wbe) « 2.46 x 1(T 9 . (5.10) 

Physically, this result may be interpreted as coming from the fact that heavy and light modes 
are interchanging energy at a constant rate, therefore rendering only a change in the overall 



amplitude of the spectrum. However, as manifest from the effective theory (4.28) the speed of 
sound is modified as: 

^ s = e^=(l + 4 V l^y\ (5.11) 

This implies the generation of non-Gaussianity noticeable in the bispectrum, as studied in 
Ref. [521E3]. 

5.2 Single turn in the trajectory 

As a next step, we consider the presence of a single turn in the inflationary trajectory. To 
simplify our analysis, we consider the specific case in which the trajectory is initially geodesic (a 
straight path), then goes through a short period in which it suffers a turn, and finally goes back 
to a geodesic state. Figure [4] shows a prototype example of such a situation. We also assume 
that throughout this process all the slow-roll parameters except for r/± remain nearly constant. 
To model this situation, we take r/± to be an analytical function of the e-fold number N: 

where AN is the number of e-folds during which the bending happens, and iVo is the e-fold value 
at which the bending is at its peak, in which case i]±(Nq) = ryj_ ma x- We recall that iV may be 
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suitably denned from conformal time r through the relation dN = ail dr. For the other slow-roll 
parameters, we choose the reference values e = 0.022 and 7711 = 0.034. These values correspond 
to a spectral index n-ji = 0.98, and to a tensor to scalar ratio r = 0.35, which are marginally 
compatible with current CMB tests [5j. Additionally, these values imply H = 10~ 5 Mpi. Figure[5] 
shows the power spectra for eight cases with different choices of the parameters AiV, r/j_ max and 
M 2 . The plotg^l contain both the spectrum obtained by solving the full coupled system of 
equations (solid line) and the spectrum obtained by solving the effective single field equation 
of motion (dashed line). For simplicity, we normalise our results in units of 2.46 x 10 -9 and 
give the scale k in units of Mpc -1 . As a reference, we have included the case 7]± = 0, which 
corresponds to the power spectrum that would be obtained in the single field case. 

The main characteristic shown by the plots are oscillatory features appearing in the spectrum. 
It may be noticed that the e-fold width AiV during which the turn takes place actually set 
the scale k of the oscillatory features. On the other hand, the amplitude of the oscillations 
is roughly dictated by the ratio Ar]± maiX H 2 / M 2 . More precisely, the amplitude of the largest 



oscillatory feature is of order SVnfPn ~ ^V±ma,xH 2 /M 2 , which agrees with the result of (5.9). 
Additionally, the match between the curve predicted by the effective theory and the full set of 
equations becomes better as M 2 /H 2 acquires larger values, irrespective of how large is /3. In 
fact, in all of the examples shown we have /3 ~ 1. 

The appearance of oscillatory features, not just a single bump, in the spectrum reflects the 
fact that both modes v and v T backreact at sub- horizon scales as the turn happens. Once 
both modes cross the horizon, the amplitude of the adiabatic mode becomes frozen (therefore 
capturing the moment in which the mode was receiving or releasing energy) while the amplitude 
of the heavy mode quickly decays due to the accelerated expansion. In fact, we have checked 
that the levels of isocurvature perturbations at the end of inflation are negligible. 

5.3 A specific example 

As a last step towards understanding the effects of curved trajectories, we discuss our results 
applied to a specific toy model, where turns are produced due to the non-trivial evolution of the 
sigma model metric. Let us consider a two-field model with fields 4> = X an d c/) 2 = tp with a 
kinetic term containing the following sigma model metric: 



lab = \ > ( 5 - 13 ) 




where T(x) is only a function of the x field and restricted to satisfy T 2 (x) < 1- The non- 
vanishing connections are T xx = —TT X /(1 — T 2 ) and T xx = T x /{1 — V 2 ) with T x = d x T, and 



3 Please note that the value of AJV in these plots corrects a factor of 2 error appearing in the published version. 
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Figure 5: The primordial power spectrum Vn{k) normalised in units of 2.46 x 10~ 9 , obtained for six different 
choices of AiV, ?7_L m ax and M 2 . The plots show a comparison between the power spectrum obtained using the full 
system of equations (solid line) and the one obtained using the effective theory (dashed line). We have chosen as 
a pivot scale the value fc* = 0.002Mpc _1 . 
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the equations of motion for the background fields are found to be 



it i r 

x - y^x 2 + 3Hx + Y^f2 V x - i^f2 V ^ =° > ( 5 - 14 ) 

^ + r^r2 x 2 + 3H ^ + i~r^ v * ~ i^y* v x =0 ' (5 - 15) 

where V x = d x V and = d^V. For concreteness, let us consider the following separable scalar 
field potential: 

%^) = %) + ^V. (5.16) 

In the particular case of T = 0, the dynamics of the two fields decouple and inflation may be 
achieved with x by a suitable choice of the potential Vo(x). If, however, T(x) 1S allowed to be 
non-vanishing for certain values of x, then a mixing between the two modes is inevitable, and 
the inflationary trajectory will be curved. Following the discussion at the beginning of Section|4j 



we choose the tangential and normal vectors T a and N a as in (4.3) and (4.4): 



-=i(x,V>), (5.17) 

(Tin V / 



N ° = <wrrpA-' p - T:i - i + r V- (5 ' 18) 

where (po = x 2 + ty 2 + ^X^- Recall that with this convention rj± is allowed to change its sign. 
The relevant background parameters describing this situation are then: 

x 2 + ip 2 + 2rxV> 

€ = 2Mf^ ' (5 - 19) 

q | xv x + ipV^ 

r?H=3+ — ^ v - — , (5.20) 

h ( x 2 + V> 2 + 2rx^ 



flVI - r 2 ( x 2 + V> 2 + 2r^ 



where -ff is given by %M^H 2 = x 2 + v + 2FxV> + 21/. For concreteness, let us consider a 
parameter T(x) having the following x-dependence 

r(x) = — o — — , (5-22) 

KX) cosh 2 [2(x-Xo)/Ax]' 1 ' 

where To is the maximum value attained by r(x). On the other hand, we take a generic potential 
Vo(x) rendering values e = 0.022 and m = 0.034 for the slow-roll parameters in the absence 
of curves. For this specific configuration, we found that the background value of e(r) remains 
nearly constant at the attractor value e = 0.022 whereas the background value of Wi(t) is more 
sensitive to the turns suffered by the trajectory, having small deviations from the attractor 
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value 7711 = 0.034. Additionally, we found two relevant time scales determining the behaviour of 
background quantities t?|| and rj±: 

Tip =M~ X , (5.23) 
T x =^ = -J^ . (5.24) 



The appearance of these time scales are actually easy to understand. First, notice that T x is the 
time during which the turn takes place whereas is the oscillation period of the massive field ip. 
We find that if T x CT^, then the background dynamics is such that 4>q = (x, ip) oscillates about 
ip = 0, meaning that both r}\\ and rj± presented oscillatory features with frequency O (t^ 1 ^ . On 
the other hand, if T x ^> T^, the background field departs adiabatically from the minima of the 
potential ip = 0, and the time evolution of 7711 and rj± is dictated by the time scale T x . This latter 
case may be interpreted as a situation where the trajectory is momentarily pushed towards one 
of the walls of the potential, as the curve takes place. Figure [6] shows the background values 
of r]± and 7711 (as functions of the e-fold number N) for the case To = 0.9, M 2 = 300H 2 and 



two values of A%, namely Ax = 0.12Mpi and = 0.084Mpi 16 In the latter case, it may be 
appreciated how the time scale appears in the shape of 7?j_. 

Figure [6] also shows the power spectrum obtained for the two described cases (right panels). 
In the present examples, the features appearing in the spectrum are not as regular as those of 
Figure [5j This is mainly because in the present situation the curvilinear trajectory contains 
several turns, in order to go back to the attractor solution. Although in this specific model the 
slow-roll parameter 771 1 appears to be sensitive to the mass scale M and the curves taking place, 
it is important to notice that this is a model dependent characteristic, and that in general t/|| 
may show various types of behaviour depending on the sigma model metric and the potential. 
In general, however, the momentary time variation of 7711 due to curved trajectories does not 
spoil the slow-roll regime, and background fields tend to quickly evolve back to the attractor 
behaviour characteristic of the single field case as soon as the bending of the trajectory stops. 
In this regard, we find that the time variation of 7711 is not relevant for the appearance of features 
in the power spectrum, and that the main contribution is that coming from the derivative 
interactions due to 77j_ in the equations of motion. 



5.4 Enhancement of non-Gaussianity 

We briefly elaborate here on another potentially observable feature not already discussed. The 
previous section computed the power spectrum of the curvature perturbation for a few examples 
where the inflaton traverses sufficiently curved regions in field space. From the results, it is 
clear that features in the spectrum will result each time the trajectory traverses a bend. These 
features are produced via the kinetic interaction between the heavy isocurvature modes and 
the light curvature mode as the turns are traversed by the background field. Crucially in these 

16 Notice that these plots are corrected versions of those appearing in the published version of this article. 
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Figure 6: Left panels: the evolution of rj± (solid line) and 5 x rju (dashed line) as functions of e-fold number N 
for two set of values of parameters Ax, To and M 2 /H 2 . In the first case, Ax = 0.12A/pi and the maximum value 
of rj± is about \rj±\ w 4, whereas in the second case, Ax = 0.084Mpi and the maximum value becomes \r)± \ ~ 6.5. 
Right panels: the resulting primordial power spectrum ■p-re(fc), normalised in units of 2.46 x 10~ 9 , obtained for 
the set of parameters used in the plots of rj±. The scale k appears in units of Mpc -1 . 



examples the heavy mode remained very massive throughout (M 2 S> H 2 ), highlighting the fact 
that heavy fields may not always be disregarded (truncated) when computing the spectrum for 
adiabatic modes. 

What is important to note is that the interaction between curvature and isocurvature modes 
implies a change in the speed of sound for the curvature perturbations- as long as M 2 3> H 2 , 
then /3(t, k) is effectively /^-independent before horizon crossing and the speed of sound may be 
written as 

cl = e~e= (l + 4r,l^\ . (5.25) 
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As is well known, a model with a speed of sound significantly smaller than unity gives rise 
to a noticeable level of non-Gaussianity of equilateral type, characterised by the non-linear 
parameter [55] : 

/fi? - 4 • (5-26) 

Thus, we are led to reason that for generic models of inflation with curvilinear trajectories in a 
multi-dimensional field space, glitches in the power spectrum are accompanied by a correlated 
enhancement of non-Gaussianity of the equilateral type, provided that the turns in the inflaton 
trajectory violate the adiabatic approximation vigorously enough- a phenomenon which we have 
argued occurs at various points in field space in many realistic realisations of inflation. Thus 
although there appear to be many models where either non-trivial modulations in the power 
spectrum (e.g. features in the single field inflaton potential [51416]) or large equilateral non- 
Gaussianity (e.g. DBI inflation j69j EO]) result, it appears that in generic multi-field models 
with curved inflationary trajectories, both are present and correlated. Evidently, the effective 
quadratic action (4.31) contains the leading higher order corrections which can also result in 
non-Gaussian signatures and implies the non-linear parameter (5.26) [66]. However to fully 
describe the bispectrum associated with the curvature perturbation, we need to properly take 
into account the cubic order action including gravity. We will discuss this issue in a separate 
publication. 



6 Conclusions 

Multi-field models of inflation contain a range of physics which goes beyond that encountered 
within the single-field paradigm. In this work we have focused on the particular case where 
all of the scalar fields remain massive during inflation except for one, which slowly rolls down 
the multi-field potential. We have found that curved inflationary trajectories can generate 
significant features in the primordial spectrum of density perturbations arising from normal 
modes becoming excited and backreacting on the dynamics of the adiabatic mode. 

To achieve these results we analysed the evolution of the quantum perturbations of a general 
multi- field setup, including the presence of a non-canonical kinetic term. Our methods are 
completely general and naturally incorporate those implemented in previous works [49^ 161] . 
where stochastic Gaussian variables are used. Moreover, although the main focus of this work 
was the study of systems where there exists a hierarchy, our results may be used to study a wide 
range of situations, including where no such hierarchies are present. 

Our formalism allows us to consider time-dependent situations beyond the regime of appli- 
cability of existing methods, such as inflaton trajectories with fast, sudden turns (regardless of 
whether the sigma model metric is canonical or non-canonical) as well as more general situ- 
ations in which the masses of the heavy fields in the orthogonal direction are changing along 
the trajectory (even if they still remain much heavier than H 2 and all other scales of interest). 
Additionally, we wish to emphasise that these non-decoupling effects have their origin in the 
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non-geodesic nature of the trajectories in field space 

Our results highlight the limitations of simply truncating heavy physics when modelling single- 
field realisations of inflation and show under which circumstances high energy effects can leave an 
imprint on the power spectrum. The main reason behind these effects is the existence of kinetic 
couplings between adiabatic and non-adiabatic modes, emerging as the inflationary trajectory 
suffers a turn. As we have seen in Section 3.1 it is always possible to change basis to a canonical 
frame where such interactions are absent. In that case, the eigenvectors of the perturbation mass 
matrix quickly vary as the inflationary trajectory turns, and we are left with the alternative point 
of view by which these high energy effects appear due to a violation of the adiabatic condition 
for truncating heavy fields. In fact, if the heavy fields are sufficiently massive, we find that we 
can construct an effective field theory for the adiabatic modes encapsulating the relevant effects 
of the full multi-field dynamics. As we have seen, such effects are not mere corrections to the 
standard single field theory, but represent entirely new contributions to the quadratic action for 
perturbations. 

We find particularly noteworthy, the presence of potentially observable signatures that result 
from a reduced speed of sound for the adiabatic perturbations during sudden turns. As a 
corollary, correlated non-Gaussianity will also manifest as a result of these sudden turns although 
a full analysis studying the details of their appearance in multi-field inflation is beyond the scope 
of this work and will be addressed in a future report [71]. Nevertheless, it would appear that 
in generic multi-field models with curved inflationary trajectories, both effects are present and 
correlated, and can potentially give information about other, much heavier, fields that would 
otherwise be inaccessible to experiment. 
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A Additional details on some background quantities 



To obtain (2.23) we proceed as follows: first, by taking a total time derivative to (2.6) we obtain 

(A.l) 



1 D 2 rh a 

90 ~3H 2 (eT a + rj a )-V (j> V a , 



'0 



dt 2 



where = T a V a . Recalling that T a = 0o/^Oj the previous equation can be re-expressed as 

D 2 T a „ , ,„ iy<t> ~ ^o) v n 



dt 2 



T a v 6 v s - v^v 



4 



-N 



(A.2) 



On the other hand, taking a total time derivative to (2.13) we may obtain yet another expression 
for the second variation D 2 T a j dt 2 , given by 



2rna 



D Z T 



( yjyjo _ Vn_ 
Wo 0o 
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V N DN a 
6n dt 



(A.3) 



Equating the last two expressions and performing some straightforward algebraic manipulations, 
we finally obtain 



DN a 
dt 



1 



Hrj±T a + — - P ab V^V b , 



H V± 



(A.4) 



where we have defined the projector tensor P' 



■ab 



1' 



ab 



T a T b — N a N b along the space orthogonal 



to the subspace spanned by the unit vectors T a and N a . That is, P ab N b = and P ab T b = 0. 
To arrive at the form of the mass matrix Cljj shown in (4.10), (4.11) and (4.12), we may start 



from the explicit form deduced out of (3.3) and (3.15) for the case of two-field models: 



2a 2 H 2 e (3 - 2tj\\ + e) 







NN 



a 2 H 2 (2-e) + a 2 V 6 



a 2 H 2 {2 - e) + a 2 V NN + a 2 H 2 M^e\ 



CItn =a 2 V ( f > N + 2a 2 H 2 r]±e , 
where we have defined 

Vd 



T a T b V a V b . 



V NN =N a N b V a V b , 
Vd. N =T a N b V a V b . 



(A.5) 
(A.6) 
(A.7) 



(A.8) 
(A.9) 
(A.10) 



Additionally K = 2M.tntn = T a N b T c N d M. abc d is the Ricci scalar parametrising the geometry of 
Ai. Notice that can be rewritten in the following way: 

V H =T a V a (T b V b ) - T a (V a T b )V b 



1 DT b 



dt 



(A.ll) 
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where, to go from the second to the third line we made use of (4.1) and relation Vn = <PoHr] 



coming from the definition of r]± in (2.18). Similarly, the quantity V^jv may be manipulated in 
the following way: 

V> =T a V a (N b V b ) - T a (V a N b )V b 

=v<t>v N - 



-Vb 



(A.12) 



where again, to go from the second to the third line, we made use of (4.1). As a final step, we 
may use Vn = <f>Hr]± to deduce: 

1 d 



V<f,V N 



dt 



-H 2 V± ( 



(A.13) 



Collecting all of these terms back into (A. 5), (A. 6) and (A. 7) we finally arrive at (4.10), ( 4.11| ) 
and (4.12). Observe that we are not able to rewrite V NN = N a N b V a V b in a similar way, since 



it involves second variations away from the inflationary trajectory. This simply means that the 
quantity Vnn must be regarded as an additional parameter of the model related to the mass of 
the transverse mode with respect to the inflaton trajectory. 



B Commutation relations for quantum multi-fields 



In this appendix we show that the commutation relations in (3.28) are fully consistent with 
the evolution of the v^(k,T) dictated by the set of equations of motion (3.31). To begin with, 
observe that in order to satisfy the commutation relation (3.28) the Af mode solutions v a {k,r) 
must satisfy the following conditions: 
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(B.l) 
(B.2) 

(B.3) 



dr dr dr dr 

To show that these relations are satisfied at any given time t we proceed as follows: first, let us 
define the following matrices: 
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I j 
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: Et^ 
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I Vvj Vv J a * 

^ dr dr 
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(B.4) 
(B.5) 
(B.6) 
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These tensors satisfy the following properties 



A IJ 


=A IJ * = 


-A JI , 


(B.7) 


B IJ 




-B JI , 


(B.8) 


E IJ 


jpij* 




(B.9) 



In other words, they are real, with A IJ and B IJ antisymmetric (E IJ has no specific symmetries). 
Because of these properties A IJ and B IJ consist of Af(Af — l)/2 independent real components 
each, whereas E IJ consists of Af 2 independent real components. Thus, in order to fix the values 
of all of these tensors we need to specify 2Af 2 — Af independent quantities. These tensors also 
satisfy the following equations of motion: 
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(B.10) 
(B.H) 
(B.12) 



0. 



(B.13) 



Taking the trace to the last equation, we obtain that the trace E = E 1 j satisfies 

dE 

and therefore E is a constant of motion of the system. Furthermore, observe that the config- 
uration E IJ = E5 IJ /Af and A IJ = B IJ = for which conditions fiO ) to pjb are satisfied 
corresponds to a fixed point of the set of equations (B.10) to (B.12). That is, they automatically 
satisfy: 

v a ij = R b ij = v eIJ = 

dr dr dr 

Therefore, it only remains to verify whether there exist sufficient independent degrees of freedom 
in order to satisfy the initial conditions E IJ = E5 IJ /Af and A IJ = B IJ = at a given initial 
time Tj. As a matter of fact, we have exactly the right number of degrees of freedom. As we have 
already noticed there exists Af independent solutions v^(k, r) to the equations of motion. To fix 
each solution v a {k,r) we therefore need to specify 2J\f 2 independent quantities, corresponding 
to the addition ofAf 2 components v^(n) and Af 2 momenta T>v^/dT(ri). However we must notice 
that the overall phase of each solution v a {k, r) plays no roll in setting the initial values for A , 
B IJ and E IJ . We therefore have precisely 2Af 2 — Af free parameters to set E IJ = E5 IJ /Af and 
A IJ = B IJ = 0. Of course, the value of the trace E is part of this freedom, and we are free to 
fix it in such a way that E/Af = 1. 

To summarise, it is always possible to choose the initial conditions for v^(k, r) and Dv^/ dr{k, r) 
in such a way that conditions (B.l) to (B.3) are satisfied. These conditions ensure the com- 
mutation relation (3.28). To finish this discussion, recall that one possible choice for the initial 
conditions for the perturbations allowing (B.l) to (B.2) to be satisfied, are precisely those ex- 
pressed in (3.34) with suitable choices for the coefficients v a (k) and TT a (k): 

v a (k)iT* a (k) - v* a (k)7r a (k) = i 



(B.15) 
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for a = 1, • • - A/*. We should emphasise however that this is not the unique choice for initial 
conditions, and in general, any choice for which E IJ = E5 IJ /N and A IJ = B IJ = will do just 
fine. 



C Effective theory for the adiabatic mode 



In this appendix we offer another deduction of the effective theory shown in Section 4.2 We 
begin by writing the action (3.16) for the particular case of two fields: 

„T\ 2 
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^) 2 -(v^) 2 -(^-C 2 ) (v N ) 2 



drd d x 



<h<r.cv N {n TN -^-2c4-) v T . 
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(C.l) 



Given that £Inn 3> |^tt| an d &nn 3* I^taH the field v N is the heavier of the two. Taking this 
scale as the scale of the heavy physics that we wish to integrate out, we can formally evaluate 
the functional integral for v N to obtain the one loop effective action for v T as 



S 



drd s x 



1 



^) 2 -(v^) 2 -(^t-C 2 )K) 2 



+ ^ drd 3 x J dT'd 3 x'0(T)v T (x,T)G{x,T;x',T')0(T')v T (x',T') + S CT ■ 
with O given by 



0(r) 



i Ltn -3 ^(7- 

ar ar 



and 



G (x, r; x' , t 



1 



(C.2) 



(C.3) 



(C.4) 



□ + n NN - e 

The term Sct renormalises the effective action for the background inflaton field, and we have to 
demand that the parameters of this effective action that satisfy the slow roll conditions rather 
than those of the bare action [72], which we presume to be the case here. In general, evaluating 
the full effective action is a highly non-trivial task. However in Fourier space, one can formally 
make the expansion: 



G{r,T',k) 



1 



1 



0. 



i + - 



where 



k +Q 



NN 



(C.5) 
(C.6) 



35 



Where implicit in the above is that if the scale M tends to the cutoff of the theory (so that 
Vnn ~ M 2 ) we can neglect the temporal derivatives in the expansion above relative to the 
mass term and the spatial derivatives (which always become significant at horizon crossing), 



thus reducing the Green's function to leading order to only the contact term ° Integrating the 



second term in (C.2) by parts results in: 



,5' 



drd d k 



' dv 
~dW 



T\ 2 



[k 2 + n(r,k)] {v T y 



with 



-i 



fi(r, k) = n 



4a 4 #V[(l + e 



i 3 H 3 rjj_(l + e- T)\\] 



to- 



(C.7) 

(C.8) 
(C.9) 



Qq is given by (4.30). Making the field redefinition ip = e^l 2 v T (and upon integrating by parts 
the resulting friction term), after some manipulations one then obtains the effective action 



S 



drd d k 



(C.IO) 



where Q(r,k ) is defined as in (4.29). Thus we see that expression (4.31) follows 
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